Gravitational Waves and 7-Ray Bursts 

Christopher S. Kochanek\°^ 

and 

Tsvi Piran^^^ 

Harvard- Smithsonian Center for Astropliysics 
60 Garden Street 
Cambridge, MA 02138, USA 

m ■ 

0> ■ Coalescing binaries in distant galaxies are one of the most promising sources 

^ ! of gravitational waves detectable by the LIGO project. They are also a co- 
^1 pious source of neutrinos, however these neutrino pulses are far too weak to 
'"^ I be detected on earth. Several years ago Eichler et al. suggested that they 
■ are also sources of 7-ray bursts (GRBs). Recently it was found that GRBs 
^ . are likely to be cosmological in origin, and coalescing binary systems^^'^"^^] are 
! probably the most promising cosmological sources. The current estimates of 

^ I the burst and LIGO rates from a cosmologically distributed population are 
ly^ ' based on the systems observed in our galaxy, t^^"^"^] These estimates are based 
^ ■ on only three binarys so there are large statistical uncertainties in the coa- 
ly-^ . lescence rate. If we accept the cosmological/coalescing binary hypothesis for 
O ! GRBs, then we get a more accurate estimate of the rate at which binaries coa- 
[ lesce and hence of the predicted LIGO signal. The association between GRBs 
^ ' and binaries can significantly improve the performance of LIGO. The detection 
p i' of gravitational radiation from a GRB source not only confirms the coalescing 
Q . binary model, but it will also provide information on the geometry and energy 

', generation mechanism of the burst. 

^ I GROt^l has shown that GRBs are located in the outer parts of the galactic halo or 

' cosmologically. In both cases all GRB mechanisms must evolve through a pair fire ball 
phase, [^^1 which makes it difficult to distinguish between types of sources based on burst 
characteristics. The association between GRBs and LIGO sources is an unambiguous, 
I verifiable prediction of the coalescing binary model, t^^l unlike gravitational lensingt^^l and 
redshift effects which are generic to all cosmological models. 

The ability of LIGO to find coalescing binaries is strongly limited by the random noise 
levels in the detectors and the need to continuously search for weak signals. There is a 
significant gain in sensitivity if the search for signals is limited to narrow windows of time 
near bright 7-ray bursts. The polarization of the gravitational waves provides evidence 
about the emission geometry of the GRB sources, because the relative strengths of the -|- 
and X polarizations is a function of the inclination angle. If the 7-rays are preferentially 
emitted along the polar axis of the binary, then the GRB sources are stronger LIGO sources 
than the average coalescing binary because binaries emit gravitational waves more strongly 
along the polar axis than in the orbital plane. 
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The evidence for a cosmofogical origin of GRBs is the combination of isotropy on 
the sky and inhomogeneity in space. The C/Cmin distribution for bright bursts has a 
—3/2 slope characteristic of a homogeneous, EucUdean distribution. Fainter bursts show 
a sharp break from this slope, indicating that the distribution becomes inhomogeneous 
at large distances. The shape of the C/Cmin for cosmological sources is determined by 
redshift effects that change the luminosity, and the effective volume. The relative 

sharpness of the break suggests the GRBs arc good standard candles, because a wide 
intrinsic luminosity distribution would have smeared out the break. [^^1 

We consider a burst with a differential luminosity per unit energy E equal to (dL/dE) = 
LqE~°'. In a detector like BATSE the number of counts, C, during a given time interval. 
At, from redshift 2; = a; — 1 is 

where is the luminosity distance and Emin is the minimum detectable energy. If the 
comoving rate of coalescences is r{z) per unit comoving volume, then the integrated rate 
of coalescences closer than redshift z = x — 1 is 

pl+z 

R{z) = 4TTr% J x-^/^Dlr{z)dx (2) 

where Dl = thDl = 2rH{x — x^^^) for the Einstein-DeSitter cosmology we use in this 
study, rn — c/Hq is the Hubble radius, and Hq — lOO/io km s~^ Mpc~^. 

For a given spectral index a and evolution model r{z), the observed distribution of 
C/Cmin is fit to infer the maximal redshift, Zmax, from which GRO detects GRBs.'^^^ 
For a constant coalescence rate, r, and spectral index a = 1, the best fit to the C/Cmin 
distribution gives Zmax = 1-45. Varying the spectral index from to 2 changes Zmax from 
2.8 to 1. The dominant uncertainty lies in the comoving coalescence rate, r{z), since a 
very strong evolution of the event rate r{z) oc t{z)^ with cosmic time t{z) can reduce Zmax 

to as low as Zmax = 0.5. 

The observed GRB rate is ~ 800 events per year^^l which corresponds to a local merger 
rate of 3.6 x 10~^hQ^ Mpc~"^ yr"^ or 3.6 x 10~^n^^ mergers per galaxy per year for a 
galaxy density of n = 10~^/iq Mpc~"^. The current rate increases by a factor of 10 if 
we assume a strong cosmological time dependence {r{z) oc t^iz)). The overall agreement 
between estimated merger rate (based on GRBs) and the estimates of 10~^^^ yr~"^ mergers 
in the Galaxy based on galactic binary pulsars^^^'^^l supports the coalescing neutron star 
model as sources for GRBs.t^^^ The agreement is better if we assume that the event rate is 
independent of cosmological time. This suggests that we use the calculated distribution of 
GRBs to predict the distribution of gravitational wave strains he, that could be observed 
by LIGO. 

The characteristic strain produced by a coalescing binary with reduced mass ^ = 
IA/insMq and total mass M = 2.8MnsMq scaled to the masses characteristic of NS-NS 
binaries, averaged over binary inclinations, and using an optimal filterl^'^-*^] with a knee at 
frequency / = 100/c2 Hz to detect the signal is 

K = 2.3 X lO-^'hof^TsMl/^f;,'/''-^. (3) 



2 



This is not the same as the instantaneous strain at frequency /, which varies as thx^^^ / Di,^ 
because the optimal filter effectively integrates the signal over many orbital periods to 
increase the signal to noise ratio. 

In Figure 1 we show the expected rate of events stronger than a given strain he for 
three different models: a constant comoving rate of r = 10~^/i~^ Mpc~^ yr"-*^ based on the 
galactic estimates, I-'^^'-'^^] the rate using the fit to the C/Cmin distribution with a = 1 and 
a constant comoving rate for which Zmax = 1-45, and the rate using the fit to the C/Cmin 
distribution with a = 1 and r{z) oc t{z)'^ for which z^ax = 0.5. The rapidly evolving 
model may be inconsistent with the local galactic estimates. 

The signal processing for the LIGO project is predicated on using two or more detec- 
tors to eliminate non-Gaussian sources of noise. Under this assumption, the instantaneous 
noise level in the LIGO detector I^'^^' is given by /i^ = fcSnifc)/ (Ff.) where fc is the 
knee frequency, (F^) = 1/5 is the angle and polarization averaged detector beam pat- 
tern, and Snifc) is the spectral density of the noise at the knee frequency fc- This is 
only the average instantaneous noise level, and to estimate the detectability of a source 
we must take into account the fluctuations of the noise. If we are using optimal filtering, 
we must make trial correlations of the filter with the detector output approximately once 
every = 0.1/^^ seconds, or about 3 x 10^/^^ trials per year. For two such Gaussian 
detectors in which we expect 3 events per year, we must have a signal to noise ratio in 
the signal of S'^/N'^ — ln(10^/^^) to be 90% confident the detectors are not seeing noise. 
This defines the detection threshold for events with an expected rate of 3 per yeart^l, 
hs/yr = (ln[10^/^'^^])^/^/in(/c) — 4:-'^hn{fc)- This estimate is insensitive to changes to 99% 
confidence rather than 90%, or 0.3 events per year rather than 3 because of the logarithmic 
dependence. Using the estimates from [3] , the initial version of LIGO should have a noise 
level of h^/yr — 8 x lO"^^/^/^ for detecting bursts with knee frequencies near / = Ifcs 
kHz. The higher knee frequency also reduces the strength of the signal by a factor of 
1.5 because of the f~^^^ dependence in he- The advanced versions of LIGO should have 
better low frequency sensitivities and noise levels of h^/y^ — 1.3 x 10~^^/c2 for detecting 
bursts with knee frequencies near / = 100/c2 Hz. The quantum limit for detecting bursts 

is hs/yr ^ 1.6 X 10-^^f-'/\ 

This is not true if we use 7-ray bursts to trigger searches for gravitational waves. As 
we can see in Figure 1, LIGO is sensitive only to signals from the rare, bright bursts, so 
we should search only a short time interval near the bright bursts for gravitational waves. 
We expect the GRB to occur between several dynamical times (milliseconds) and several 
cooling times (seconds) after the coalescence of the binaries. If we are conservative and 
search for one minute before each of the ten brightest bursts in a given year, then we are 
conducting only 6000/^"^^ trials rather than 10^/^"^^ trials, and the detection threshold for 
7-ray burst events is 1.5 times larger than the h^/yr estimates. An increase in the detection 
threshold by a factor of 1.5 increases the rate by a factor of three. Events correlated with 
GRBs not only come at a known time, but also from a known direction. This reduces the 
number of trial correlations further, but it is not as important as restricting the search in 
time because of the poor spatial resolution of LIGO. 

If the use of two detectors does not eliminate all non-Gaussian errors in the LIGO 
signal processing, then the use of Gaussian statistics to estimate h^/yr will fail catastroph- 
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ically for the random search mode, because the noise estimate is working way out on the 
tail of the Gaussian distribution. If the noise has a non-Gaussian tail, triggered searches 
using the GRBs may be the only way LIGO can detect coalescing binaries. 

Coalescing binaries emit gravitational waves much more strongly along the polar axis 
than in the orbital plane, with oc (1 + cos^ i) and /ix oc 2 cosi where i is the inclination 
angle of the binary. I^^J The characteristic amplitude he oc + l/ixP)"*^^^ includes an 

average over the inclination angle between the observer and the binary. If the 7-ray 
emission from a coalescing binary is beamed along the polar axis, we have underestimated 
the average gravitational wave sig nal by the factor (l+lla^/ie+llxVie+x^/ie+x^/ie)^/^ 
where x — cosip and is the opening angle around the normal into which the bursts are 
beamed. If the GRB is associated with a disk formed during the mergerl^^'^^J then we 
would expect the emission to be beamed along the polar axis. The maximum gain is a 
factor of 1.6 in and the gain is 1.2 for '0 = 30°. Because the detection volume increases 
as the cube of the detection limit on the strain, the rate increases are factors of 4 and 1.9 
respectively over unbeamed bursts. If the bursts are detected by LIGO the beaming can 
be detected from the polarization of the bursts: unbeamed bursts tend to be dominated 
by the + polarization, while beamed bursts arc dominated by the x polarization of the 
gravitational waves. If the bursts are beamed along the polar axis, then the GRB rate 
misses fraction cosi/' of the coalescences. 

The hypothesis associating GRBs with coalescing binaries has several important impli- 
cations for gravitational wave searches using LIGO. The GRB rate can accurately estimate 
the rate of coalescences as a function of the gravitational wave strain on the earth. The 
LIGO sensitivity must have a strain sensitivity he greater than 10~^^"^/io for an event rate 
of one per decade, and greater than 10~^^'^/io for a rate of once per year. The strain sensi- 
tivity of LIGO is considerably greater if the brightest GRBs serve as triggers to search for 
gravitational waves in a narrow window of time near the burst, and in the direction of the 
burst. The strain sensitivity is roughly 50% greater in a coincidence experiment compared 
to a random search, and the detection rate is roughly three times higher. 

The presence or absence of coincidences between GRB and LIGO events may be the 
final proof or disproof of the coalescing binary hypothesis. Once the strain sensitivity 
reaches 10~^^ '^ho LIGO can begin to set limits on GRB models, initially by setting limits 
on the redshift distribution. If the sensitivity reaches 10~'^^-^ho then either LIGO begins 
to detect coalescing binaries in coincidence with GRBs, or the coalescing binary model 
is dead. If they are related, then polarization of the gravitational waves can be used to 
determine the emission geometry of the GRBs. The time delay between the merging of 
the binary and the GRB will help to determine the emission mechanism. The shortest 
possible delays will be comparable to the dynamical time (milliseconds), and longer delays 
might be viscous or weak interaction time scales. 

The burst rate is an underestimate of the LIGO rate. For a given class of binaries, 
beaming of the 7-rays may mean that only a fraction of the coalescences are seens as 
GRBs. Moreover, not all types of coalescing compact binaries may cause bursts. The 
galactic rates are only estimated from NS-NS binaries, and there may be a comparable 
number of NS-BH mergers and a much smaller number of BH-BH mergers. I^^^ It is not 
clear whether NS-NS or NS-BH binaries are more likely to produce bursts, although the 
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dynamical stability of orbits and evolution time scales favor NS-NS binaries over NS-BH 
binaries J^^^ 
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Figure 1: The figure shows the integral event rate of binary coalescences with character- 
istic strains greater than he- The three models are a fit to the GRB rates with a constant 
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comoving rate and a = 1 spectral index (solid line), a fit to the GRB rates with a comov- 
ing rate proportional to r{z) oc t{z)'^ (dashed line), and a constant event rate of 10~^h~^ 

Mpc~'^ yr~^ (solid with points). The strain scales as ^oA*]vs-^Af5/c2^^^- T^^^ estimated 
LIGO sensitivity hs/yr is shown by the vertical solid line, and the sensitivity improvement 
from using coincidences with GRBs is shown by the vertical dashed line, strains he is 
shown as a function 
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